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Abstract 

We consider a bivariate rational generating function 

under the assumption that the complex algebraic curve V on which Q vanishes is smooth. For- 
mulae for the asymptotics of the coefficients {Ors} are derived in [PW02] . These formulae are 
in terms of algebraic and topological invariants of V, but up to now these invariants could be 
computed only under a minimality hypothesis, namely that the dominant saddle must lie on 
the boundary of the domain of convergence. In the present paper, we give an effective method 
for computing the topological invariants, and hence the asymptotics of {ors}, without the min- 
imality assumption. This leads to a theoretically rigorous algorithm, whose implementation is 
in progress at |http : / /www . mathemagix . org| 
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1 Introduction 

Consider a power series F(z) = ^^.arz"', where r varies over integer vectors in the orthant (Z+)'^ 
and z"" denotes the monomial ■ ■ ■ . We say that F is the (ordinary) generating function for 
the array {or}- In analytic combinatorics, our aim is to derive estimates for Or, given a simple 
description of F as an analytic function. An apparatus for doing this for various classes of function 
is developed in a series of papers |PW02[ IPW041 IBP08[ IPemlO] : see also the survey |PW08] . Most 
known results on multivariate asymptotics concern either rational functions or quasi-powers. In the 
case of rational functions F(z) = P{z) /Q{z), the analysis centers on geometric properties of the 
pole variety V := {z : Q{z) — 0}. 

When V has singularities, tools are required such as iterated residues, resolution of singularities 
or generalized Fourier transforms, and analyses exist only for specific examples. When V is smooth, 
if certain degeneracies are avoided, asymptotic formulae for Or may be given in terms of certain 
algebraic and topological invariants of V. Denote normalized vectors by 



Asymptotic formulae for Or depend on the direction r in which r is going to infinity. For example. 
Theorem 3.9 of |PemlO[ Theorem 3.9] gives an asymptotic formula for Or in terms of a sum over a 
set S of quantities that are easily computed via standard saddle point techniques. The set S = S(r) 
is a subset of the set saddles of saddle points of the function h : V ^ M. defined by 



i=i 

The set saddles is readily computed but membership in the subset S is not easily determined, and 
in fact there is no known algorithm for doing so. 

The main result in this paper is a characterization of 5 when d — 2 and the complex algebraic 
curve V is smooth. This leads to a completely effective and rigorous algorithm for asymptotically 
computing Or. Precise statements of these results are somewhat technical, so we refer to the main 
text. The organization of the paper is as follows. 

SectionOreviews background results from elsewhere which reduce the computation of asymptotic 
formulae to identification of the set S and of some path segments through each c G S. In particular, 
Section [3 begins with the integral representation of the general coefficient a^, Section [TT] defines 
the residue of a meromorphic form. Section 12.21 reduces the integral for to a lower-dimensional 
integral with some parameters yet to be specified, and Section 12.31 selects a chain of integration for 
this integral that results in an explicit asymptotic estimate in terms of saddles, S and the so-called 
critical height c* (Theorem l2.13p . 
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The new material begins in Section [31 Beginning with the topology of V near the coordinate 
axes (Theorem 13. 2p . we give a topological characterization of the critical minimax height c* to 
which the cycle of integration can be lowered (Theorem I3.4[) and of the set S of contributing saddles 
(Theorem 133]). Section S] demonstrates how these topological computations of c* and S can be made 
completely effective fAlgorithm 14.81 and Theorem 14. 9p . Details of implementation are given, as well 
as a discussion of uniformity and boundary cases. 



2 Background 

In this section we review the framework for deriving estimates for a^, beginning with results for 
general F, then specializing to rational functions, and smooth pole variety V. This section may 
be skipped by readers who want only to understand the effective computation and not 
the underlying analysis. Most of the material in this section is valid for any number of variables. 
When we need to assume a bivariate function, we will state this but will also change the notation 
to use X and y in place of zi and Z2 and (r, s) in place of r. 

Computation of via complex analytic methods begins with the multivariate Cauchy integral 
formula. 

(27rz)'* Jrp zi ■ . . . ■ Zd 

Here the torus T is any product of a circles in each coordinate sufficiently small so that the product 
of the corresponding disks lies completely with the domain of holomorphy V of F. This version of 
Cauchy's formula can be found in most textbooks presenting complex analysis in a multivariable 
setting, and follows easily as an iterated form of the single variable formula; see, for example, |Sha921 
page 19]. The integrand is may be written as z~'^uj where 

F(z) 

uj:^ujf-= — — dz. (2.2) 

Zl- Zd 

Let (C*)'' (C\0)'' denote the set of complex d-vectors with all nonzero coordinates and denote 
V' := V n (C*)''. A key step in obtaining asymptotic estimates for is to transform the Cauchy 
integral equation (27ri)'^ar = /j, z^'^up via an identity valid for any meromorphic form with a simple 
pole: 

/ z-^uj = 2TTi / Res (z"'"a;) (2.3) 

Jt J a 

where Res (•) is the residue operator defined in Section [2. II and a is the intersection class defined in 
Section [2?2] Putting these together and specializing to d — 2 leads to Lemma [2^ below: 



^ /x-V^Res(.), 
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2.1 The residue form 



The residue form Res {rj) is a holomorphic (d— l)-form on V'. The specification of this form will not 
be important for the subsequent analysis, but for completeness we include the following definition. 
If 77 = (P/Q) dz is any meromorphic form on a domain U with simple pole on a set V C [/, then we 
define 

Res{r]):=L*e (2.4) 

where 9 is any solution to 

dQA9 = Pdz. (2.5) 

Existence and uniqueness are well known, and a proof can be found in [DeVlOi Proposition 2.6] 
among other places (see, e.g., |AY83| ). One important property of the residue form is that if 
g : V — C is holomorphic then 

Res {g-n)^ g- Res (77) . 
In particular. Res {z^^ uj) — z^'"Res (w). 



2.2 The intersection class 

Because Res (z^"" lj) is holomorphic on V', the integral over a. {d— l)-cycle i? C V' depends only on 
the homology class of B in Hd-iiV'). The chain of integration a in (|2.3|) is really a homology class 
known as the intersection class of T with V'. The construction of the intersection class is quite 
general and may be found in the literature. Conditions are given in |PW121 Appendix A] for the 
intersection class to be uniquely defined. When it is not, however, the integral over a must be the 
same for any choice of intersection class, a, so we will not pursue it further here. 

Because we will need an explicit construction of a cycle in this class, we will give a quick con- 
struction of a and derivation of (12. 3p . We begin with a form of the Cauchy-Leray residue theorem, 
which may be found in [DeVlOl Theorem 2.8]. 

Lemma 2.1 (Cauchy-Leray Residue Theorem). Let rj he a meromorphic d-form on domain U 
with pole variety V C [/ along which r/ has only simple poles. Let N be a d-chain in U , locally the 
product of a [d — I) -chain C on V with a circle 7 in the normal slice to V, oriented as the boundary 
of a disk oriented positively with respect to the complex structure of the normal slice. Then 

/ ri — 27ri / Res (77) . 
Jn Jc 

□ 

Corollary 2.2. Let F,P,Q and V he as above and let rj be meromorphic. Let Aii and M2 be 
compact d-manifolds in and let H be a smooth homotopy from Mi to Mi 2 that intersects V 
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transversely in a compact manifold C. Then 

/ V ^ V — 27ri / Res (77) . 
JMi Jm2 Jc 



□ 



Proof: Let N be the boundary of a tubular neighborhood W olC in C^. Then 77 is holomorphic on 
H\W . The boundary of H\W is A^2 ^ A^i and dr] vanishes on H\W because the differential of 
any holomorphic d-form vanishes on C*. Tubular neighborhoods of a smooth algebraic hypersurface 
always have a product structure (see, e.g., |PW12[ Proposition A.4.1]). Therefore, Stokes' Theorem 
and Lemma 12.11 together yield 



/ 1] = I V + V ^ V + 27ri / Res (ry) . 
JMi Jm2 Jn JM2 Jc 



□ 



Remarks 2.3. One sees that the orientation of C must be chosen so that the boundary of its product 
with a positively oriented disk in the normal slice is homologous to Mi — M.2- Also note that result 
is true when H\W \s any cobordism: we used only that d{H \ W) = — M.i + N , not that H 
was a homotopy. 

Let Te be the torus {|zi| = • • • = \zj\ = e} and let T^x be the torus {\zi\ ^ • • • = = 
e , \zd\ — L}. Fix e > and suppose that for sufSciently large L, the torus does not intersect V. 
We claim that for r^; sufficiently large, 

/ z-^oj^O. (2.6) 

Indeed, 



|z"''w| < £-^.=1 '•^L-'-'' sup |w(z)| 

where f3 is the degree of the rational form u). The chain of integration has volume 0{L), from which 
we conclude that Jj, ^ tT^ u) — 0(l?-^^~'^'^\ Once r^j > 1 + /?, we therefore have J^, ^ uj = o(l) 



as i — cx). But Te^i avoids V for sufficiently large L, hence by Stokes' Theorem, Jrp ""a; i 



IS 



constant for sufficiently large L, hence zero, establishing p.6p . 

The intersection class may now be constructed as follows. Let be the homotopy defined by 
HL{z,t) = (zi, . . . , Zd-i, (1 - t)zd + t{L/e- l)zd) . 
Suppose that H^^^] intersects V transversely. We let 

a = [C] where C := iJL[T,] n V . (2.7) 
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This is independent of L once L is sufficiently large. We take the orientation of C as in the remark 
following Corollary 12.21 

Summing up the residue and intersection cycle constructions we have: 
Lemma 2.4 (Residue representation of Or). Suppose d — 2. Then for sufficiently small e > 0, 

"^=(2^?^ l^x-'^y-^ResH (2.8) 
holds with Res defined as in (|2.4p - (|2.5p and the intersection class a defined in (|2.7p . 

Proof: Let S be the set of x e C such that there is a sequence {xn,yn) in V with and 
yn — >■ oo. The set S is finite and therefore min{|x| : a; € S', x 7^ 0} is strictly positive. Taking e 
less than this minimum guarantees that for sufficiently large L, the torus T^^l does not intersect 
V. The set of {x,y) G V where dQ/dy vanishes is also finite. Again, taking e less than the least 
nonzero modulus of such a point guarantees that the homotopy intersects V transversely. Having 
verified these two suppositions, the construction of a is completed. Now apply Corollary 12.21 with 
Mi=T^, M2 = T^.L, N — C and 77 = z^'^uj, where e is small enough so that C T). Using ()2.6|) to 
see that 77 = and combining with Cauchy's integral formula gives 

{2TTifar = / z"'"a; = 2ni / Res {z^^uj) = 2ni / z""" Res (w) . 

□ 

Remark 2.5. The only use thus far of the restriction to d = 2 is in verifying that Hl intersects V 
transversely and T^,l not at all. Trans versality may be accomplished for any d via a perturbation 
and it is not hard to handle a nonempty intersection of V with T^^l, so no essential use has been 
made yet of the restriction to d = 2. 



2.3 Morse theory 

Saddle point integration techniques require that we choose a chain of integration C in the class a to 
(roughly) minimize the modulus of the integrand z^'^uj. It will suffice to minimize the factor of z^'', 
the form oj being bounded on compact chains. The log modulus of z^"" is log \zj\. Letting 

= kl = X]j=i ''j; have 

d 

log I z""" I =^rj\og\zj\ ^nht{z) 

which explains our definition of h in (|l.ip . It will be helpful to consider the manifold V from a 
Morse theoretic viewpoint with height function h. 
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Definitions 2.6. Let JH be any manifold and let /i : — > M be any smooth, proper function 
with isolated critical points (a critical point is where Vft. vanishes). Let V-"^ denote the subset of 
V consisting of points z with h{z) > c. Define V'"^, V-'^, V^"^ and V^"^ similarly. Define the set 
saddles = sadd\es{M,h) to be the set of critical points for h. The set ^.[saddles] is called the set of 
critical values. For fixed (r, s), the set saddles is a zero-dimensional variety, whose defining equations 
we call the critical point equations: 

Qix,y) - (2.9) 
sxQx-ryQy = 0. 

Here and henceforth, we use subscript notation for partial derivatives. 

When : V — > K is the height function /if above, by convention we extend h to V taking 
h{0, y) — h{x, 0) — +0O, so that points of V \ V are in every V-"^. 

Theorem 2.7. Let ft, : — > K. 6e a smooth, proper function with isolated critical points. 

(i) Suppose there are no critical values of h in the interval [a,b]. Then the gradient flow induces 
a homotopy in M between and A^-°, and also between M-'^ and Ai-"^. This homotopy 

carries M=^ = dM^'' = d^^'' to M=°- . 

(ii) Let a < c < b such that c is the unique critical value of h in [a, b]. Let cri, . . . , cr^ enumerate the 
critical points of h at height c. Then H^{M.-^ . M.-"") decomposes naturally as the direct sum 
of H^,[Ai^^ , M^"") , where M^j is the intersection of A4 with a sufficiently small neighborhood 
ofaj. 

(Hi) Let C be a cycle in A4 and let c* be the infimum of c such that C is homologous in to a 
cycle in Then c is either — oo or a critical value, and in the latter case, for sufficiently 

small e > 0, the cycle C projects to a nonzero homology class in H^,{M,M-'^~'^). 

Proof: Part {i) is the first Morse lemma, a proof of which may be found for example in |Mil63l 
Theorem 3.1] or |GM88| . A generic perturbation separates the critical values, one may then prove 
part (m) by a standard excision argument. Part {Hi) is proved in |PemlO[ Lemma 8]. □ 

Remark 2.8. Classical Morse theory combines the first Morse lemma (no topological change between 
critical values) with a description of the change at a critical value (an attachment map). The results 
we have quoted require only the first Morse lemma. For this it is not necessary to assume h is a 
Morse function, but only that it is a smooth, proper map with isolated critical points. 

Example 2.9 (bivariate harmonic case). Suppose d — 2. Let cr be a critical point of /i on V at 
height c and suppose that locally h is the real part of an analytic function. Then there is an integer 
k > 2 and there are local coordinates ip : U ^ V on a neighborhood of the origin in such that 
-0(0) = cr and 

h{^{z)) = c+ Re{z''}. 
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Figure [T] illustrates this for fc = 4. The shaded regions are where h > c. On the right, the connected 
local components of {/i > c} are number counterclockwise from 1 to and the local components of 
{h < c} are numbered correspondingly. The pair Ai := (V-'^+'^ H ip[U],V-'^^'^ D il^[U]) is homotopy 
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Figure 1: A saddle of order 4 and the path 71^2 

equivalent to a disk modulo k boundary arcs. Let denote an arc from the boundary to the center 
that remains in the j*'' local component of {h < c}. Then Hi{A4) is generated by the relative 
1-cycles 7i — jj. We have the relations J2i j ^i-jili ~ Ij) = if and only 

j J 

for every i < fc. In particular, {7^ — 7^+1 : l<i<fc — l}isa basis for Hi{AA). 

We remark that our assumptions imply that h has no local extrema on Ai. Therefore, as c 
decreases, the connected components of M.-'^ can merge but new components never arise. 

This leads to a canonical form in dimension d = 2 for classes in iJi(V"), choosing a representing 
cycle so as to be in position for saddle point integration. 

Proposition 2.10. Let h = hf. Letj be any cycle in iJi(V'). Then either^ is homologous to a cycle 
supported at arbitrarily small height or there is a critical value c, and a nonempty set ^ C saddles 
of critical points at height c, such that the following holds. 

The cycle 7 is homologous in HiiV') to a cycle C in V-*^*; the maximum of h on V is 
achieved exactly on S; in a neighborhood of each cr £ S 0/ order k — k{a) > 2, there is 
a subset paths(<T) of {1, . . . ,k{a) — 1} such that the restriction of C to the neighborhood 
is equal to J2ieT{a) 7* - It+i ■ 

As promised, integrals of the residue form Kes {x^'''y^^ujp) are known. To summarize this, we 
let F — P/Q he rational, with uj defined in (j2.2[) and Res (w) defined in Section I^TTl With ji as 
in Figure [1] we then have the following explicit integration estimates. Note that the magnitude of 
x~'^y~'^ is exp[(r + s)c*] when (a;,y) G S, which makes the remainder term exponentially smaller 
than the summands. 
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Lemma 2.11 (saddle point integrals). Let a — (a;o, j/o) G saddles, satisfy the critical point 
equations (|2.9p . Denote the order of vanishing of h at (a;o,?/o) by k — fc(a;o,yo)- Then 
(27ri) Xy._^.^j Kes {x~'^y~'^u!) is given by the formula $((T, i) where 

$(a, .f^^''''^ ^ x,^y,^ + O (exp[(c. - e)(r + .)]) (2.10) 

^27r(r + s) 

where K{a, i) is a constant term given by p. lip when d = 2 and more generally by an explicit formula 
which may he written in a number of ways, e.g. WW02\ Theorem 3.5] or \BBBP08\ Theorem 3.3]. 
When k = 2, the remainder term is uniform as {xQ^y^^r, s) varies with (xo,2/o) G saddles(r, s) and 
{xQ,yo) varying over a compact set on which k is constant. When k > 2, the estimate is uniform 
for fixed (xo,?/o) o-nd {r, s) od with r = Xs + 0{1), where X = xQx/{yQy) evaluated at (xo,?/o)- D 

Remarks 2.12. 

(i) In fact a uniform estimate holds for r = As + 0(s^/^) although this is more complicated and 
involves Airy functions. 

(ii) When d = 2 the constant K{a,i) is given in (PW021 Theorem 3.1] as 

K{a,^) = P{x,y)J^-±^^ (2.11) 
y s xip 

evaluated at (xo,2/o): where ip{x,y) is the expression in equation (j4.4p below. 

(iii) When fc > 2, the constant K(a, i) is the leading term in a series expansion for A near a. The 
explicit formula involves partial derivatives of P and Q up to order k and may be derived from, 
e.g., |PW021 Theorem 3.3]. 

Putting together Lemma with Proposition l2 . 101 and Lemma r2 . 1 1 1 leads to the following asymp- 
totic expression for Or, which has appeared several times before; for example, when a is a minimal 
point, this goes back to |PW021 Theorem 3.1]; for an alternative formulation see [BBEPOSl Theo- 
rem 3.3]. 

Theorem 2.13 (smooth point asymptotics). 

^ J2 ^<y,^)+0{exp[{c,-e)ir + s)]) . 
iepaths(<T) 

□ 

This completes the summary of known results. We now turn to the main work which is to identify 
and compute the quantities c«,S and paths (cr) appearing in Theorem l2.13l and formula (12.101) . 
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3 Identifying and S 



We begin by constructing the intersection class. We will take a homotopy which keeps the torus 
small in the x-coordinate. If V intersects the y-axis transversely at (0, y), then the set V n {|a;| < e} 
will have a component that is a topological disk nearly coinciding with a component of V-'^ with 
c = — f log e. To see that something like this holds more generally, we recall some facts about Puiseux 
expansions. Consider the multivalued function x i— >■ yix) that solves Q{x,y(x)) — 0. When a; is in a 
sufficiently small punctured neighborhood of 0, there is a finite collection of convergent series of the 
form 

yix) = J2 ^1^'^' (3-1) 

j>jo 

where jo G Z, A; € and Cj„ =/= 0; such a series yields k different values of y for each x; together, 
the collection of series yields each solution to Q{x,y) — exactly once. We always assume that k 
is chosen as small as possible. A proof of these facts may be found in [FS09[ Theorem VII. 7]. Our 
assumption that (5(0,0) 7^ implies jo < in every Puiseux series solution. Those with jo — 
correspond to points at which V intersects the x-axis, the value being equal to the coefficient cq. 
Under the assumption that V is smooth, for each such y- value there will be only one seriefjfl. If fc = 1, 
the intersection with the y-axis is transverse but it is possible that k > 2 and V is smooth but with 
non-transverse intersection. Finally, series with jo < correspond to components of V n {|a;| < e} 
with y going to infinity as x goes to zero. We begin by showing that each series gives a single 
component and describing the embedding of the component. First, to clarify, we give an example in 
which all these possibilities occur. 

Example 3.1. Let Q{x,y) — 1 — 3y + 2y^ — 6x1/ + x^y^. It is easy to check that V is smooth. 
The total number of solutions should be the y-degree of Q, namely 5. Setting x — we find two 
points (0, 1) and (0, 1/2) of V on the y-axis. These correspond to two series solutions with jo = 
and fc = 1: 

y{x) = 1 + 6a; + 72x'^ H 

1 3 45 2 
= 2-8"+32" +••• 

There are also two series solutions with jo < 0, one with fc = 1 and one with fc = 2: 

/ \ 1-1/2 3 19v3 1/2 

y{x) = - + • • • . 

The former of these yields two solutions. Although we do not need it here, we remark that these 
five solutions can be enumerated by counting y-displacement along the Northwest boundary of the 



®In fact the analyses in this paper are valid when self-intersections of V are allowed at points (0, j/), with the result 
that the closures of the components of V n {|x| < e} may intersect. 
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Newton polygon: going up from (0,0) to (0,2) counts the first two solutions, going from (0,2) 
to (1,4) yields the two solutions y ~ ±l/-\/3x and going from (1,4) to (3,5) yields the solution 
y ~ 6x~^. 

Let X-" (respectively F-'^ denote the union of those components of V-" containing arbitrarily 
small values of x (respectively y). Say that the direction of a Puiseux series solution to Q{x, y{x)) = 
is —jo/k. The direction is always nonnegative because jo is nonpositive due to (5(0,0) ^ 0. Let 

bad := {/3 > : /3 is a direction of a series solution to Q{x, y{x)) = 0} 

denote the set of positive directions. 

Theorem 3.2 (V' near the axes). Fix r with r/s ^ bad. Suppose that h has isolated critical points. 
Then when e> is sufficiently small and c is sufficiently large, the following hold. 

(i) The graph of each Puiseux expansion of y{x) over the punctured disk {0 < |a;| < e} is a 
component ofV (1 {\x\ < e}. 

(a) Each such component is a topological disk and its boundary winds k times around the y-axis. 

(Hi) The components of X-'^ are exactly the graphs over some neighborhood of of those Puiseux 
series whose directions are less than r/s. 

(iv) Components ofV' fl {|a;| < e} and X-" with the same Puiseux series are homotopic in V'. 

Proof: Fix any Puiseux series solution X^j>j(j CjX^^^ on a punctured disk U and define g{x) := 
^Cjcc-' . From this we define the function 

X ^ {x^,g{x)). 

We will show that G is a diffeomorphism onto its image. For sufficiently small e the images of the 
punctured e-disk under different Puiseux series are disjoint, which will finish the proof of (i) and {ii). 

First we check that G is one-to-one on a sufficiently small punctured disk. If not then we may 

find arbitrarily small xi and X2 such that x\ = X2 and g{xi) — g{x2). This means that there is some 
jj,h jQQ^ Qf unity ^ and such that g{x) = g{^x) for arbitrarily small x. For some m, the functions 
x"^g{x) and x"^g{^x) are holomorphic, and if they agree on a sequence converging to zero then they 
must define the same function on U. It follows that their Taylor expansions agree, and hence that 
the only nonzero coefficients are those indexed by a multiple of k. This contradicts the minimality 
of k in the Puiseux series representation and we conclude that g is one-to-one on a sufficiently small 
disk. For any smooth function (f), the function {x'',(f){x)) h- > x is diffeomorphism near any .x 7^ 
because x*^ 1— >■ x is locally smooth except at 0. This proves that G is a diffeomorphism onto its 
image, establishing (i) and (ii). To continue, we need a lemma. 
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Lemma 3.3. When e is sufficiently small, for any 9 the radial arcs 

{te'\y{te'%), 0<t<€ 

are monotone in height, that is, 

he{t) :-/^(^e'^y(^e^«)) 

is monotone increasing or decreasing. Furthermore, as t ^ 0, hg{t) — > ±oo according to whether the 
direction (3 =^ —jo/k of the Puiseux series is less or greater than r/s. 

Proof: Write 

y{x) = cx~^[l + <l){x)) 

where the remainder satisfies 4'{x) — o(l) and x(j)'{x) = o(l) as x — >■ 0. Plugging this into the 
function H{x, y) :~ — f log x ~ s log y gives 

d sP-f (l)'{x) 

—H{x,y{x)) - 



dx X 1 + 4>{x) 

1 (^o ^ a;</)'(a;) 
— — sp — r 



X \ 1 + ipix) 

s/3 — f 



(3.2) 



X 

because we have assumed sf3 ~ f ^ 0. The function h is the real part of H, which along with the 
chain rule yields 

-rheit) = cos 6* Re-^H(x,y(x)) — sin 6' Im -^H(x,y(x)) . 
at dx ax 

Plugging in (I3.2p then gives 

dhg{t) sfi — f 
~~dt 1^ 

which finishes the proof of the lemma. □ 

Proof of Theorem 13.21 continued: Each component B of X-^ has, by definition, points with 
arbitrarily small values of x which arc therefore on the graph of one of the Puiseux series. Taking 
a sufficiently small neighborhood, the graphs of the series are disjoint, therefore the map from 
component to series is well defined. We need to see that the range of the map from components to 
series is exactly those with directions les than r/s. If B has direction less than r/s then ft, — > oo 
as a; —> on the corresponding component, whence the corresponding component of V' H < e} 
maps to B. On the other hand, if the direction of B is greater than r/s then h — oo as a; — > 0, 
from which it follows that X-'^ cannot contain points in B with arbitrarily small x values and hence 
such a component is not the graph of a Puiseux series about a; = 0. This proves {Hi). 

Let B be a component of X-'^ where c is sufficiently large so that the projection onto x is 
contained in a sufficiently small disk so that Lemma 13.31 applies. A radial homotopy in the x- 
coordinate (one that maps x := pe^^ to f{9)x and maps y{x) to y{f{9)x)) then shrinks B to the same 
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component with any specified greater value of c, or the corresponding component of V H {|a:| < e} 
with as long as e is at most the in-radius of the original x-neighborhood. This establishes (iv) and 
also that either type of component remains in the same homotopy class as c or e is varied. □ 

Switching the roles of x and y in Theorem I3.2[ it follows that for sufficiently large c, the set 
V-^ is the disjoint union of X-"^ and Y-'^. Define Cxy to the least value of c for which this remains 
true, which is also that greatest value of c such that V-'^ has a component containing both a point 
(0, y) and a point {x, 0). By the first Morse lemma, if h is smooth and proper, the topology of V-"^ 
cannot change between critical values of h. In particular Cxy must be a critical value Cj. The next 
two results complete the description of and S. 

Theorem 3.4 (characterization of c*). Assume that V is smooth and that h = hf is smooth and 
proper with isolated critical points. Then 

(i) For each c > Cxy, the cycle dX-'^ is in the homology class a from Lemma \2.4\ 

(ii) C-^ — Cxy ■ 

For each critical point a with h{a) = c* let k = k{a) denote its order, let A/" denote a small 
neighborhood of it in V. Label the components oi JV O {h > c,} in counterclockwise order by 
Ri, . . . ,Rk as in Figure [TJ Each such region R will be called an a;-region if i? C X^"^ and a y-region 
if i? C y^*^. Because c > Cxy, these are mutually exclusive. 

Theorem 3.5 (characterization of S and paths). 

(^jj CT G S if and only if the regions Ri, . . . , Rk include at least one x-region and at least one 
y-region. 

(ii) Let Y{j) — 1 if Rj is a y-region and Y(j) = otherwise. Let Ca denote the relative cycle in 
iJi(A/'n V''*+%A/'n V*"'^) represented by 

k 

Then represents the projection of a to this homology group, from which we easily recover 
paths(cr). 

(Hi) The cycle C* := X^tre-'^o' represents the projection of a to Hi{V'^*^'^ . 

□ 

Proofs: Let ci > • • • > c^, = c^y be the critical values of h that are at least Cxy. Let W denote the 
set of c G M for which dX^"^ e a. To establish (i) we will show: 
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(a) If c G W and q > c > Ci+i then the entire closed interval [q+i, is a subset of W; 

(b) If c e W and c> ci then [ci, oo) C W; 

(c) If Ci e W and i < i' then q — e G W for all sufficiently small e > 0; 

(d) The set W contains some c > ci . 

By part (i) of Theorem 12.71 the spaces : cl > t > q+i are all homotopic in V', hence 

represent the same homology class in iJi(V'). By continuity, this is true as well for t — c-i and 
t — Ci+i, which proves (a). The proof of (b) is identical. 

To prove (c), we first assume that there is a single critical point a with h(a) = Ci. Let k >2 he 
the order to which the derivatives of h vanish at a. Recall the neighborhood U from Example 12.91 
and the analytic parametrization ip -.U -^V with h{il){z)) ~ Ci + Rejz''}. Locally, the image of U 
is divided into 2k sectors with h > Ci and h < Ci in alternating sectors. Figure [5] shows V-'^ (shaded) 
for three values of c in the case k ~ 2. A circle is drawn to indicate a region of parametrization 
for which h = Ci + Refz*^}. In the top diagram c > Ci, in the middle diagram c — Ci and in the 
bottom, c < Ci. The arrows show the orientation of 9V-'^ inherited from the complex structure of 
V. The pictures for k > 2 are similar but with more alternations. Consider the first picture where 




Figure 2: V- and its boundary for three values of c 

c > Ci. Because c > Cxy, each of the k shaded regions is in X-'^ or Y-'^ but not both. Let us term 
these regions "a;-regions" or "y- regions" accordingly. Because Ci > Cxy, this persists in the limit as 
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c -I Ci, which means that either all k regions are in Y~ or all k regions are in In the former 

case, dX-'^ does not contain tr for c in an interval around Ci and the first Morse lemma shows that 
ig^>ci+e Nomotopic to dX-'^^~'^. In the latter case we consider the cycle dX-'^^^'^ + dB where 
B is the polygon showed in Figure |31 Because we added a boundary, this is homologous to 
But also it is homotopic to dX'^^'^'^: within the parametrized neighborhood, the lines may be shifted 
so as to coincide with dX'^^^'^ , while outside this neighborhood the downward gradient flow provides 
a homotopy. In the case where there is more than one critical point at height we may add the 
boundary of a polygon separately near each critical x-point, that is, each point that is a limit as 
e 4, of points in 



Figure 3: V-^*+' and V-^' ^ differ locally by a boundary 



Finally, we recall from Lemma [2.41 that a cycle C representing a is constructed as Hl\T^] n V. 
Here e is sufficiently small and L sufficiently large that this set is precisely {(a;,y) € V : |a;| = e}. By 
part [iv) of Theorem 13. 2[ the union of these homotopic to dX-"^ for sufficiently large c. Checking 
that the orientation given by Remark l2.3l is the same as taking boundaries of regions oriented by the 
complex structure of V, we conclude (d). Part [i) of Theorem l3.4l is a direct consequence of (a) - (d). 

Part [i) implies that for any c > Cxy, there is a cycle in the class a that is supported on V-*^, 
which shows that c* < Cxy. To see that c* > Cxy, we suppose not and argue by contradiction. The 
construction in [PemlOl Lemma 3.8] shows that for each Ci > c,, there is a cycle representing a 
supported on V-'^'"'"'^ and any such cycle projects to zero in the relative group (V,V'^*~'^). For a 
contradiction, it suffices to check that dX'^''^'^ does not project to zero in i?i(V, V-'^*"'). 
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Suppose first that there is a single critical point a at height Cxy. Recall from Example 12.91 that 
H^^iy^^y^"^ ^\)-^^y~^') is generated by {7^ — 7i+il < i < fc — 1}. The only relation among these is that 
the sum of all of them vanishes. It follows that a sum of some subset of these vanishes if and only 
the subset is all or none. Referring back to Figure |31 it is clear that the cycle dX^^'^y is represented 
in Hiiy^'^y^^ , v<cxy-e') |;j^g gyjj^ gf |;j^g ^.^^^ _ Qygj- thosc 1 loT wHch the shaded sector between 
these two path segments is an x- region. By the definition of Cxy, the point u is a limit point of both 
x-regions and y-regions. In other words, the subset of shaded sectors that are in X^'^'^y is a proper 
subset of all fc shaded sectors. It follows that dX^^'^y does not vanish in _ffi(V, V-'^'=*'~'^), and we 
have our contradiction. This establishes (m) of Theorem l3.4l in the case of only one critical point at 
height c* and also gives the characterization of Ca in part (ii) of Theorem 13.51 in this case. 

When there is more than one critical point a with h(a') — c^, it follows from the definition of 
Cxy that for at least one such cr, the fc local regions include at least one x-region and at least one 
y-region. The formula in part (m) of Theorem 13.51 continues to hold at each a. This vanishes when 
all the regions are of one type, giving the characterization of 5 in part (i) of Theorem 13.51 The 
direct sum decomposition in part (m) of Theorem 12.71 allows us to write C* as X^o-e" ^cr, establishing 
part (iii). □ 



4 Effective computation 

By describing c^,5 and in terms of x- and y-regions. Theorems 13.41 and 13.51 give us a means to 
compute c*,S and C*. One possible approach is as follows. Instead of following 3V-'^ down from 
c — oo to c*, we follow paths upward from each critical point until we find what we're looking 
for. A critical point of order fc has k computable steepest ascent directions. The backbone of our 
computation will therefore be: 

1. Compute the critical points and order them by decreasing height. 

2. For the highest saddle, for each of the k ascent paths, check whether it converges to the x-axis 
or the j/-axis. 

3. if there is at least one of each type, set c* equal to the height of this saddle, set Co- according 
to part (m) of Theorem 13.51 and do the same for any other saddles at this height. 

4. If all are of the same type, then continue to the next lower saddle and iterate. 

It should be clear that no new theorems are required in order to implement this program. However, a 
number of computational apparati are required in order to make such a program completely effective; 
this section, which contains no major theorems, is devoted to implementation. 
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Chief among the needed computational apparati is a way to compute ascent paths. Discretizing, 
we have the problem: given {x,y) S V, produce a point {x' ,y') with h{x' ,y') > h{x,y) such that 
there is a strictly ascending path in V from {x,y) to {x',y'). Assuming we can do this, we must 
also ensure that this process ends after finitely many steps with identification as an x-component 
or a y-component. This, together with some computer algebra, will complete our preparation for 
implementation. 

4.1 Ball arithmetic and Mathemagix 

For our implementation we have chosen the platform Mathemagix [vdHLM+02] . This platform 
was designed for rigorous computations with objects of both algebraic and analytic nature. On the 
one hand, the algorithms in this paper indeed rely on purely symbolic pre-computations, which will 
be done using Grobner bases. On the other hand, our ascent paths are a special case of rigorous 
computations with analytic functions. Interval arithmetic is a systematic device for carrying out 
such computations |Moo661 IAH83[ INeuQOl IMKCOQi IRumlOj . We will use a variant, called ball 
arithmetic [vdHOQj . which is more suitable for computations with complex numbers. The use of 
interval arithmetic to achieve full rigor contrasts to other homotopy continuation methods such as 
the package Bertini |SW07) which use bootstrap testing that is extremely reliable but still heuristic. 
Although there is a large community for interval arithmetic, these techniques are not so common in 
other areas. For this reason, we will recall the basic principles of ball arithmetic. 

Ball arithmetic 

Given c £ C and r E M-°, denote by B{c,r) the closed ball with center c and radius r. In what 
follows, a complex ball number is a ball B{c, r) with c E Q[i] and r E Q-*^. We think of such a ball 
number as representing a generic complex number z E B{c,r). The standard arithmetic operations 
can be defined on ball numbers as follows: 



Given * E {+, — , •} and B{c",r") — B{c,t) * B{c',r'), these definitions have the property that 
z E B{c, r) and z' E B{c', r') => z * z' E B{c", r") . 



Similar definitions can be given for other operations, such as division, exponentiation, logarithm, 
and so forth. Notice that any element z E Q[i] gives rise to a ball number B{z, 0). 

In practice, Q is usually replaced by Fp.g, the set of floating point numbers x = m2'^ whose 
mantissa m S Z and exponent e G Z have bounded bit lengths p eN and q eN respectively. In that 



Bic,r)+Bic',r') 
B(c,r)~B(c',r') 
B{c, r)-B{c',r') 



B{c + c',r + r') 
B{c-c',r + r') 
B{c-c',i\c\+ry +r\c'\) 
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case the operations in Fp.^ have additional rounding errors, and the definitions of the operations on 
ball numbers must be adjusted to take into account these additional errors. Furthermore, we must 
allow for the case where r = -\-oo. For more details, we refer to |vdH091 Section 3.2]. 

Ball arithmetic allows for the reliable evaluation of more complicated expressions (or programs) 
/(zi, . . . , Zk) that are built up from the basic operations. Indeed, by induction over the size of the 
expression, it is easily verified that 

zi e B{ci,ri), . . . , Zfe e S(cfe, Tk) /(zi, . . . , Zfe) € /(S(ci, ri), . . . , i3(cfc, Vk)) ■ 

Furthermore, it can be shown that ball arithmetic remains continuous in the sense that the radius 
of f{B{ci, ri), . . . , B{ck, rk)) will tend to zero if ri, . . . , all tend to zero. 

Complex algebraic numbers 

Thus far, we have no reliable zero test for ball numbers. In one direction, if B{c, r) n B{c', r') = 0, 
then we can be sure that z ^ z' for any z S S(c, r) and z' G B{c' , r'). However the converse does not 
hold: if B{c, r) n B{c' , r') ^ 0, then we do not know whether z — z' ioi z € B{c, r) and z € B{c' , r'). 
Of course, for some algorithms it is necessary to be able to decide equality. Fortunately, this can be 
accomplished for the subclass of algebraic numbers. 

Represent a complex algebraic number z by a triple {P, c, r) where P e Q[t] is a square free poly- 
nomial and B{c, r) is a ball number such that z is the unique root of P in S(c, r). This representation 
is of course not unique. The condition that z be the unique root of P in B{c, r) may be replaced by 
a more explicit sufficient condition due to Krawczyk and Rump |Kra69[ IRumSOi IRumlO) . Consider 
the expression 

B{c', r') = $(S(c, r)):=c-^^ + [l- ^^^^^) • B{Q, r), (4.1) 

evaluated using ball arithmetic. If B(c',r') is contained in the interior of B{c,r) then it is certified 
that P admits one and only one root z € B{c, r). In particular, z = if and only if P is divisible 
by t. Moreover, for any e > 0, there exists an iterate B{cn,rn) — <i?"(S(c, r)) such that r„ < e and 
z € B{cn,Tn)- In othcr words, the triple {P,c,r) may be replaced by a new triple {P,Cn,rn) for 
which the radius r„ is arbitrarily small. 

We will denote by A the set of complex algebraic numbers. Using the above representation, 
arithmetic operations in A are easy to implement. For instance, assume that we want to compute the 
sum of two algebraic numbers represented by (P, c, r) and (Q, c', r'). We first use symbolic algebra to 
compute the square free part R of the annihilator of all sums of a root of P and a root of Q. Replacing 
B{c, r) and B{c' , r') by smaller balls if necessary, we next ensure that B{x", r") B{c, r) + B{c' , r') 
satisfies the above Krawczyk-Rump test for R. The triple {R,c",r") then represents the sum. 
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Combining Grobner basis techniques, e.g., |CL098| . with algorithms for complex root finding of 
univariate polynomials jSch82[ IBFOO] . the following can be shown in a similar way: there exists an 
algorithm which takes as input a zero-dimensional ideal of Q[ti, . . . and produces as output all 
complex solutions in the form of vectors of ball numbers describing disjoint poly-balls each containing 
a unique root. 



Ascent paths 



Given a point zq G with f{zo) 7^ 7^ f'{zo), it will be important to be able to compute a line 
segment of ascent, namely another rational complex number zi such that |/| increases on the line 
segment from zq to zi. First, we note the following procedure to check whether a ball lies entirely 
in the sector 5^/4 of complex numbers whose arguments are strictly between — 7r/4 and 7r/4. We 
consider a ball B{c,r) that is known to intersect and we let B{c',r') denote Re {B{c,r)^} in 
ball arithmetic. If r' < |c'| then zero is not in B{c',r'), hence by continuity S(c, r)^ lies in the open 
right half-plane. Because B{c,r) is known to intersect 8-^/4, it follows that B{c,r) lies in rather 
than In other words, to conclude that B{c, r) C Sj^/^, it is sufhcient that r' < |c'| in the ball 

number that is the real part of the square of B{c, r). 

Proposition 4.1. Let f be locally analytic such that f and /' can be evaluated using ball arithmetic. 
Let Zq be a rational point at which f{zo) 7^ 7^ /'(^o)- 



(i) For e > sufficiently small, the function 
evaluated at u = S(0, e) is contained in 8-^/4. 



(4.2) 



(ii) For such an e > 0, the function \f\ is strictly increasing on the line segment from zq to 

Remark 4.2. Let eo{f,ZQ) denote the supremum of those e > satisfying (j4.2p when evaluated in 
ball arithmetic. Although implementation-dependent, this quantity is continuous as a function of 
Zq. It is also easy to find an e > which satisfies (|4.2|) and which is larger than eo{f, zq)/2: just 
keep doubling e as long (|4.2|) remains satisfied. 

Proof: Change variables to the function g defined by 



1 f/^ ^ 



Thus 5(0) = (?'(0) = 1. Because 1 is in the interior of iS'^/4, continuity of ball arithmetic implies that 
g'{B{0,€)) lies in 6*^/4 for sufficiently small positive e. The function in (|4.2I) is precisely g' , which 
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proves (i). For (ii), observe that for u G [0, e] we have 

9{u) = / g'iv)dv e S^/i. (4.3) 



This imphes Aig {g'{u)/g{u)) G (— 7r/2,7r/2) for all u G [0, e]. Thus \g\ increases strictly on [0, e], 
and changing variables back to / gives (ii). □ 

More generally, we will need to compute a segment of ascent for algebraic points where the 
derivative vanishes to some order k. The same argument as for Proposition 14.11 easilv shows the 
following. 

Proposition 4.3. Let f be locally analytic such that f and its first k derivatives may be evaluated 
using ball arithmetic. Let zq be a rational point at which /(zq) =/= /'-'^■'(•^o) while /'-"'^'(^o) = ■ ■ ■ = 
= 0. Define 

1 

where any choice of the fc*'* root is allowed. Then when e > is sufficiently small, ^'■'^^(^(O, e)) C 
St^^4 and for such an e, the magnitude of f will increase strictly on the line segment from zq to 
zq + e[f{zQ)/ f^''^ (zq)]^/'' . Furthermore, computing such an e for each choice ofl/k power and taking 
the minimum ensures that \ f\ increases on all k line segments simultaneously. □ 





r/(^o)i 


i/fc\ 


\zq+U 


[/'(^o)J 





4.2 The use of ascent paths to compute invariants on V 

We now harness these computational devices to compute the topological invariants of V that are 
required for asymptotics, namely c*,S and paths. We begin by transferring ascent segments of a 
function (Propositions 14. l l and 14. 3p to ascending paths on the Riemann surface V. Let W denote the 
set of points z = (a;, y) G V \ saddles with x G Q[i] and dQ/dx{z) ^ 0. We may represent elements 
of W as ball numbers with first coordinate radius zero and second coordinate radius small enough 
that the ball contains only one root of Q{x, •). 

Proposition 4.4 (rigorous ascent step). There is a ball- computable function (j) : W ^ W with 
the following properties. Let Zq = (a^Oiyo) G W and denote zi :— {xi,yi) := (/'(zq). Then the line 
segment [xq, xi] lifts uniquely to a curve in V connecting Zq to zi, along which h is strictly increasing. 
Furthermore, <j) may be chosen so that /i(zi) — ft.(zo) is bounded below by a positive constant ck on 
any bounded subset K of W whose closure avoids saddles and points where dQ/dx vanishes. 

Proof: Given Zq G W, let y be the locally analytic function such that y{xo) = yo and Q{x, y{x)) = 0. 
Apply Proposition 14.11 to /(z) := exp[h{z,y{z))] and Zq = Xq, obtaining a segment [a;o,a;i] = 
[zo,Zo + ei/(zo)//'(^o)- It is easy, computationally, to choose ei always to be at least eo/2. Define 
4>i{zq) :— zi :— {xi,yi). The lifting of [a;o,a;i] to V is a path along which h increases, hence 
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h{(j>i{z)) — h{z) is strictly positive. Let us check that this difference is bounded below by a positive 
constant on compact sets. 



We know ei is bounded away from zero on compact sets avoided by saddles and dQ/dx because 
eo is. Letting h{zo,t) denote the derivative of |/| along the lifted segment at zq + tf{zo)/ f (zq), and 
observing that h{zo,0) is bounded below by a positive constant on compact sets avoiding saddles, 
we conclude that /i(0i(z)) — h{z) is indeed bounded below by a positive constant on compact 
sets, K. 

This gives the conclusion we desire except that xi may not be rational. To correct this, note 
that all we used about ei was that was guaranteed to be in the interval [ceo,eo] for some constant 
c. Having computed ei, we may choose e > 0.99ei such that ef(zQ)/f'(zo) is rational. Denoting 
the endpoint of the new shorter lifted path by (/'(z), the same argument as before now shows that 
h{4>{z)) — h{z) is bounded from below by a positive constant, finishing the proof. □ 

Having defined a step from a point not near a saddle, we next find a way to ascend out of a 
saddle. The following proposition describes a way to do this. The proof is nearly identical to the 
proof of Proposition 14.41 and is omitted. 

Proposition 4.5 (ascent from a saddle). Let zq = (a;o,yo) G saddles with xq G Q[«]- Suppose 
dQ/dx does not vanish at zq. Then we may compute a rational ')jj{xo) such that the union of the 
radial line segments 



The following proposition allows us to terminate an ascent path when it comes close enough to 
the a;-axis or y-axis. 

Proposition 4.6 (arrival at X or Y). Fix f with r/s ^ bad and let h = h}. Let Cmax denote the 
greatest critical value of h and Cmin denote the least critical value. Let e he small enough so that there 
are no critical points {x,y) with \x\ < e. Then any point {x,y) £ V with \x\ < e and h{x,y) > c,nin 
is in X'"^"""' . 

Remark. The hypothesis is formulated so as to be easily checked using ball arithmetic: letting J be 
the ideal generated by the critical point equations (12. 9p , we take e to be sufficiently small such that 
the ball solutions to J are disjoint from the ball B{0, e). Also note that the roles of x and y may be 
switched to yield an analogous result. 

Proof: Let e be as in the hypothesis. We know that the graph of V over ,8(0, e) is a union of 
graphs of Puiseux series. Recall that as x — on one of these components, ft — oo except possibly 
for some "low" components on which y is unbounded and h — >■ — oo. By hypothesis, h cannot take 
on a critical value on any of these components, so we have h > Cmax on the high components and 




lifts uniquely to a union of paths from Zq on V on each of which h is strictly increasing. 



□ 
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h < Cmin on the low components. This accounts for all solutions to |x| < r on V', hence \x\ < r 
together with c > Cmin suffices to assure that a point (x, y) G V is in X^'^"'"'. □ 

The next proposition allows ascent paths to terminate when they come sufficiently near a saddle. 
Without this improvement there is a danger that an ascent path could converge to a saddle along an 
infinite sequence of ever smaller steps. This necessary improvement is also a big time-saver. Step 2 
of the algorithm at the beginning of Section 2] calls for us to identify which axis is the limit of each 
ascent path from a given saddle. Proposition 14.71 classifies an ascent path once it comes near any 
higher saddle, thereby saving the remainder of the journey to the axes. 

Proposition 4.7 (bypassing a saddle). Let c > be a critical value of h and let (xo,?/o) be a 
critical point at height c in X-'^. Suppose there is an e > for which the following conditions hold. 

(i) For each x such that \x — X{)\ < e there is at most one solution to Q{x, y) = with \y — yo\ < e; 
(ii) {xo,yo) is the unique critical point {x,y) of h onV with |a; — xo|, |y — j/o| l£ £• 

Then \x — xq\, \y — yo\ < e and Q{x, y) = together imply (x, y) £ X^'^' . 

Remarks. An identical result holds if X-'^ and X^'^' are replaced by Y-'^ and Y^^' throughout. 
Also, hypotheses (i) and [ii) are easily checked in ball arithmetic. 

Proof: Our first hypothesis on e implies that the set 

S■.^Vf^{\x~ xo\ < e} n {\y - yo| < e} 

is a graph over B{xQ^e) of a univalent holomorphic function. Our second hypothesis implies that 
V^^' does not intersect S. Because V^^* is the boundary of V^"^' this implies that the entire set S 
is contained in one component of V^"^*. Because (a;o,yo) is an element of S and is also in X^"^*, we 
conclude that S C which is the desired conclusion. □ 

Let (Ti , (72 , • ■ • , cTm bc ball algebraic representations of the critical points in weakly descending 
order of height. Fix f ^ bad. The algorithm outlined at the beginning of Section |4] may now be 
described in detail, proved to terminate and proved to produce a correct answer. In the end, the 
index jo will be such that Cj„ — c«; each critical point will be classified as x, y, mixed or unclassified; 
the set S will be the set of mixed saddles; the unclassified saddles will be the ones lower than c*. 

Algorithm 4.8 (Computation of c*,S and C*). 

1. (initialization) 

• Set j := 1. 

• Let e be as in Proposition |4]6] and let e' be the analogous quantity with x and y switched. 
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• For 1 < j < m set class(crj) := "unclassified". 

• Set c* := — oo. 

2. (main loop) Repeat while h{aj) > c*: 

• Set Zq :— aj. 

• Apply Proposition 14.51 with {xo,yo) — zq for each of the k choices of power and let 
Zi for 1 < i < k denote the other endpoint of the resulting line segment. 

• For z = 1 to fc, initialize subclass^ to "unclassified". 

• For i = 1 to fc, initialize done to False and repeat until done: 

— If the x-component of Zi has modulus less than e then set subclass; :="x" and done := 
True; 

— else if the y-component of Zi has modulus less than e' then set subclass^ :="?/" and 
done := True; 

— else if Zi is in the ball B{ap) in the conclusion of Proposition 14.71 for any {xo,yo) = 
ap : p < j then set subclass^ SubclasSp and done := True; 

— else set Zi :— (j){zi) as in Proposition 14.41 

• — If clasSi ="a;" for all 1 < i < fc then set class(crj) ="a;". 

— else if clasSi ="?;" for alll <i < k then set class((Tj) ="?/". 

— else set class((jj) :— "mixed". 

• If c» = "undefined" and class((Tj) = "mixed" then set c* := h{aj). 

• 3 j + 1- 

3. (finishing the computation) 

• Set S = 0. 

• For i = 1 to j, if class((Ti) = "mixed" then do: 

— S:=SU{a,}. 

— Compute via Theorem 13.51 with Y{i) ~ ^subclass -y' 

Theorem 4.9. This algorithm terminates and correctly computes c*,S and {Co- : a G S}. 

Proof: The sequence Zq, Zi,(j){zi), (f>{(j){zi)) , . . . defines a polygonal path in x whose lifting to V is 
an ascent path for h along which the value of h on each step increases by an amount bounded away 
from zero. We must therefore reach the terminal conditions of the done loop in a number of steps 
bounded by L, where L = M — Cmin divided by this minimum step increase and M is the greatest 
height of any {x,y) e V with |a;| = e. The k{(j) ascent paths from each critical point a verify 
whether each of the k regions is and cc-region or a j/-region; the algorithm computes Cxy and thus, 
by Theorem 13. 4| computes c» . Theorem 13.51 shows that the remainder of the algorithm computes S 
and C* . □ 
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4.3 Uniform asymptotics as f varies 



Let 

monkey := {r : has a degenerate critical point } 

where degenerate means having order 3 or higher. We use this name because the degenerate saddle 
of order 3 is called a "monkey saddle" (downward regions for both legs and the tail) . 

Proposition 4.10. Assume for nondegeneracy that Q is not a binomial. Let J be the ideal in 
C[x, y, A] generated by the following three polynomials. 

Q 

xQx ~ ^yQy 

yQliQx + xQxx) + xQliQy + yQyy) - IxyQxQyQxy (4.4) 

Then {r,.s) S monkey implies r/s is the X-coordinate of a solution to J. In particular, the set 
monkey is the intersection o/(0,oo) with the zero set of the elimination polynomial f of X in J. 



Proof: The first two generating polynomials for J define an algebraic function 

A(„,:.g (4.5) 

on the variety {Q — 0}. This is the direction function of |PW02j in the sense that a point (x, y) E V 
is critical for /if. if and only if X{x,y) — r/s. The derivative of A on V with respect to, say, x, is 
given by 

which vanishes exactly when the Wronskian XxQy — XyQx vanishes. Using the explicit formula (j4.5p 
for A shows the Wronskian, up to some factors of y and Qj,, to be the third of the given generators of 
J. Vanishing of the derivative of A on V is a necessary condition for the coalescing of two solutions 
to A = c on V. This proves the main conclusion. If there is no elimination polynomial for / in J 
then A is constant on V which implies that Q is binomial. We have assumed not, from which the 
last statement follows. □ 

Let A :— bad U monkey denote the set of troublesome directions. We see from Proposition 14. 101 
and the definition of bad that A is readily computed. Its complement in the arc of f in the positive 
quandrant is a finite union of intervals with endpoints whose slopes are algebraic. On each interval 
in A'^, the asymptotic estimate for Or in Lemma 12.111 is uniform over compact sub-intervals. This 
allows us effectively to give uniform asymptotics in all nondegenerate directions. 

Algorithm 4.11. 
• Compute A 
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• Enumerate the intervals of A'^ in the open arc parametrized by A G (0, cxi) and choose one 
rational value r{I) in each interval /. 

• for each interval / do: 

— Set f = r{I) 

— Use Algorithm 14.81 to compute c* , S and C* 

— Use Lemma 12.111 to write a uniform estimate for Or with f e / by summing (|2.10p over 
cr e S and i S paths((T). 

• Interpret the formula as holding uniformly over compact subintervals of /, where multi-valued 
quantities depending on f are extended by homotopy from their values at r{I). 

Discussion of remaining directions 

To complete the asymptotic analysis of ar, we need to know what happens in the remaining direc- 
tions. For each f G monkey, we have a formula given by inserting the estimates from Lemma [2.111 
into Algorithm 14.81 These estimates are known to extend to Airy functions on the rescaled window 
r = As + 0(s^/^). Because the focus of this paper is how to make the computation of S and paths 
effective, we do not discuss the extensions of p.lOp to Airy function estimates in the present work. 

When r/s = Aq G bad, there is a projective point of V of finite height. It is possible that this 
point is a smooth point and not a critical point, in which case a change of chart maps gets rid of 
it; the formulae for the two adjacent intervals / and /' will agree and will be valid throughout the 
union / U /' U {Ao}. More often, however, the projective point fails to be a smooth point. In this 
case the formulae on / and /' will in general be different and asymptotics in the direction Ao will be 
given by something other than (|2.10p . 
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